Saturation effects in pp scattering in the impact-parameter 
representation* 

O.V. Sclyugin 1 ^ and J.R. Cudell* 

The impact of unitarity is considered in different approaches to saturation in impact-parameter space. 
The energy and momentum-transfer dependence of the total and differential cross sections and of the ratio 
of the real to imaginary parts of the scattering amplitude are obtained in a model that includes soft and 
hard pomeron contributions, coupled to hadrons via the electromagnetic form factor, ft is shown that the 
hard pomeron may significantly contribute to soft physics at the LHC. A similar conclusion can also be 
reached in the framework of non-linear approaches to unitarisation of the BFKL pomeron. 



1. Introduction 

One expects that non-linear effects will 
enter the BFKL equation in the non- 
perturbative infrared region, i.e. at large im- 
pact parameters. This is a different regime 
from that connected with the Black Disk 
Limit (BDL), in which saturation occurs at 
small impact parameters first. 

The common viewpoint is that saturation 
will lead to a decrease of the growth of a to t- 
But the estimates of the energy after which 
saturation will be important vary a lot be- 
tween different models. 

Unitarity of the scattering matrix is con- 
nected with the properties of the scattering 
amplitude in the impact parameter represen- 
tation, which is equivalent at high energy to a 
decomposition in partial waves. The scatter- 
ing amplitude can then saturate the unitarity 
bound for impact parameters b < bi. To sat- 
isfy the unitarity condition, there are different 
models. Two of them are based on the solu- 
tion of the unitarity equation SS^ = 1. First 
of all, in the [/—matrix approach pQ, one ob- 
tains a ratio <J e il&tot — * 1, as s — > oo. The 
second possible solution of the unitarity con- 
dition, which is the usual one, corresponds to 
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the eikonal representation 

/>oo 

T(s,t) = i bdbJ (bA)[l-exp(-x(s,b))] (1) 
Jo 

with t = A 2 . If one takes the eikonal phase 
in a factorised form %(s, b) — h(s) fib), one 
usually supposes that, despite the fact that 
the energy dependence of h{s) can be a power 
h(s) ~ s A , the total cross section will satisfy 
the Froissart bound a to t < C \og\s). 

We find in fact that the energy dependence 
of the imaginary part of the amplitude and 
hence of the total cross section depends on the 
form of f(b), i.e. on the s and t dependence of 
the slope of the elastic scattering amplitude. 

1.1. Eikonal representation and atot 

Let us first take a Gaussian for the eikonal 
phase /(&) ~ exp(— b 2 /R 2 ). As was made first 
by Landau, let us introduce the new variable 
y = exp(— b 2 /R 2 ). We can then calculate 
the integral exactly and obtain that 

T{s : t)~iR 2 [r(0,s A ) + A logs] , (2) 

where T(a,z) — t a ~ x e _< dt and, in our 
case, F(0, s A ) — > as s — > oo . If R 2 is 
independent from s, we have atot ~ l°g( s )7 
whereas if R 2 grows not faster then log(s), 
the total cross section becomes proportional 
to the Froissart bound a to t ~ l°g 2 ( s )- 

However, let us now take a polynomial form 
f(b) ~ s A /b 4 . Such a form comes, for exam- 
ple, from the BFKL equation. In that case, 
we can introduce the new variable y = 1/6 4 
and obtain 

f°° 1 

T(s,t) ~ i / — [1 -exp(-s A y)] 

Jo vVv 

= 2v^s A/2 . (3) 



2 



So, in this case, the scattering amplitude even- 
tually violates the Froissart bound! 

If we introduce an additional small constant 
radius r which removes the singular point b = 
in f(b) and take f{b) ~ s A /[b 4 + r% the 
answer, after some complicated algebra, is 



2r 4 



2r 4 



•(4) 



The asymptotic value of the Modified Bessel 
functions is 

J 0)1 (, A /(2r 4 )) ~ r 2 / (^ A/2 ) • (5) 

We again obtain for asymptotically high en- 
ergies 



T(s,t) ~ i V^s A/ 7 (2V2) 



(6) 
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Figure 1. The total cross section of proton- 
proton scattering calculated in the eikonal 
representation (hard line: with an exponen- 
tial form; dashed line: with a Gaussian form) 
compared with a log 2 (s) dependence (dash- 
dotted line). 



Finally, let us consider the more compli- 
cated (but most interesting) case, of an ex- 
ponential form for f(b) ~ exp(— mb). The 
corresponding amplitude in the t representa- 
tion is 

T(s,t) ~ i s A q F p [(l, 1,1), (2,2,2), -s A ]. (7) 

with p = q = 3, and where the function q F p is 
the hyper-geometric function. 

A numerical estimate of this integral shows 
that we obtain for the exponential form of the 



eikonal an energy dependence of the scatter- 
ing amplitude in the ^-representation which 
can be approximated as T(s,i = 0) ~ 
i a log 2 (s/so), with large coefficients a — 4.5 
and s = 135 GeV 2 . 

Hence, in this case, the scattering ampli- 
tude obeys the Froissart bound, but with a 
large scale so and a large coefficient a. This 
leads to a weak energy dependence at mod- 
erate energies followed by a fast growth at 
super-high energies. 

Of course the eikonal representation, which 
most of the time leads to a unitarity an- 
swer at every impact parameter, leads to am- 
plitudes which do not saturate at finite en- 
ergy: the eikonal representation for the scat- 
tering amplitude in 6-space, in the form 1 — 
exp(— x(s, &)), reaches the BDL only asymp- 
totically. However, this representation is not 
the only possibility, and it may be more useful 
to consider the effects of saturation by consid- 
ering parametrisations in s and i, transform- 
ing them to impact parameter space, and im- 
posing directly the BDL as an upper bound 
on the amplitude in s and b. 

1.2. Non-linear effects and o t ot 

A different approach to saturation is found 
in the studies of the non-linear saturation pro- 
cesses, which have been considered in a per- 
turbative QCD context |2I3| . Such processes 
lead to an infinite set of coupled evolution 
equations in energy for the correlation func- 
tions of multiple Wilson lines 0]. In the ap- 
proximation where the correlation functions 
for more than two Wilson lines factorise, the 
problem reduces to the non-linear Balitsky- 
Kovchegov (BK) equation |4I5| . 

It is unclear how to extend these results 
to the non-perturbative region, but one will 
probably obtain a similar equation. In fact we 
found simple differential equations that repro- 
duce either the U-matrix or the eikonal rep- 
resentation. We can first consider saturation 
equations of the general form 



dN(Z,b)/d£ = S(N) 



(8) 



with TV the true (saturated) imaginary part 
of the amplitude. We shall impose that (a) 
TV — > 1 as s — > 00, (b) dN/dt; — > as s — > 00, 
(c) S(N) has a Taylor expansion in TV, with 
the hard pomeron Nb are — f(b)s A as a first 
term. This enables us to fix the integration 
constant by demanding that the first term of 
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the expansion in s A reduces to Nb are . 

Inspired by the BK results, we shall use 
the evolution variable £ = logs. If we want 
to fulfil condition (c), then we need to take 
S(N) = AN + 0(N 2 ). Conditions (a) and (b) 
then give S(N) = A(N - N 2 ) as the simplest 
saturating function. The resulting equation 

dN/dlogs = A(N- N 2 ) (9) 

has the solution 

N = f(b)s A /(f(b)s A + 1) (10) 

One can in fact go further: eq. 10 has been 
written for the imaginary part of the ampli- 
tude. If we want to generalise it to a complex 
amplitude, so that it reduces to © when the 
real part vanishes, we must take: 

dA/dlogs = A(A + iA 2 ) (11) 

The solution of this is exactly the form 
obtained in the U-matrix formalism, for 
ZU( S ,b) = s A f(b). 

Many other unitarisation schemes are pos- 
sible, depending on the function S(N). Wc 
shall simply indicate here that the eikonal 
scheme can be obtained as follows: 

<9iV/<91ogs = A(l-iV)(-log(l-iV)) (12) 

Other unitarisation equations can be eas- 
ily obtained via another first-order equation. 
The idea here is that the saturation variable 
is the imaginary part of the bare amplitude. 
One can then write 

8N/dN bare = S'(N)=> 
dN/d\ogs = [dN hare /d\ogs] S'(N) (13) 

with N^are the unsaturated amplitude. 

This will trivially obey the conditions (a)- 
(c) above, and saturate at N = 1. Choos- 
ing S'(N) = 1 — N gives the eikonal so- 
lution whereas S'(N) = (1 — N) 2 leads to 
the U-matrix representation (|l(Jfl . We have 
thus shown that the most usual unitariza- 
tion schemes could be recast into differential 
equations which are reminiscent of saturation 
equations |4I5| . Such an approach can be used 
to build new unitarization schemes and may 
also shed some light on the physical processes 
underlying the saturation regime. 

2. Conclusion 

In the presence of the hard Pomeron [H], the 
saturation effects can change the behaviour 



of some features of the cross sections already 
at LHC energies. Some forms of the eikonal 
phase in the factorising eikonal representation 
can lead to a violation of the Froissart bound. 
Non-linear effects which work in the whole en- 
ergy region supply an acceptable growth of 
the total cross sections. Saturation leads to a 
relative growth of the contribution of periph- 
eral interactions. The most usual unitariza- 
tion schemes could be recast into differential 
equations which are reminiscent of saturation 
equations Such an approach can be used to 
build new unitarization schemes and may also 
shed some light on the physical processes un- 
derlying the saturation regime. 
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